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ABSTRACT 

We present a method for constraining the evolution of the galaxy luminosity-velocity (LV) relation 
in hierarchical scenarios of structure formation. The comoving number density of dark-matter halos 
with circular velocity of 200 km s -1 is predicted in favored CDM cosmologies to be nearly constant over 



the redshift range 



5. Any observed evolution in the density of bright galaxies implies in turn 



a corresponding evolution in the LV relation. We consider several possible forms of evolution for the 
zero-point of the LV relation and predict the corresponding evolution in galaxy number density. The 



Hubble Deep Field suggests a large deficit of bright (My < —19) galaxies at 1.4 



2. If taken at face 



value, this implies a dimming of the LV zero-point by roughly 2 magnitudes. Deep, wide-field, near-IR 
selected surveys will provide more secure measurements to compare with our predictions. 



Subject headings: cosmology: theory 
universe 

1. INTRODUCTION 



Within the framework of hierarchical structure forma- 
tion, galaxies form within virialized dark-matter (DM) 
halos. The halo properties and their time evolution can 
be predicted robustly within a given Cold Dark Matter 
(CDM) model using N-body simulations or analytic ap- 
proximations. A fundamental property of a DM halo is 
its maximum circular velocity, which is likely to be related 
to the observable internal velocity of the luminous galaxy 
that resides in the halo. 

In the local universe, it is well known that galaxies obey 
fundamental scaling relations, including the relation be- 
tween total luminosity and internal velocity (commonly 
known as Tully-Fisher or Faber- Jackson relations), which 
we will focus on here. This LV relationship surely reflects 
some fundamental aspect of galaxy formation, although 
its origin is poorly understood. An observational determi- 
nation of the redshift evolution of this relationship would 
provide crucial clues as to its physical origin and would 
constitute an important constraint on models of galaxy 
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formation. However, a direct observational determination 
of the internal velocities of galaxies at high redshift is ex- 
tremely difficult. Rotation curves have been obtained for 
a small number of galaxies up to redshift z ~ 1 (Vogt et 
al. 1996, 1997, 2000). These results suggest that there 
is little evolution in the Tully-Fisher relation for objects 
similar to local large spirals out to z ~ 1 (but see Mallen- 
Ornelas et al. 1999). At higher redshifts (z > 1), no direct 
observations are yet available. 

While galaxies of different morphological types obey dif- 
ferent LV relations, one can define a mean LV relation for 
any given morphological mix. This can be used to make 
a direct connection between the distribution functions of 
galaxies in terms of their internal velocity and luminos- 
ity. Clearly, if the evolution of the velocity function is 
predicted by theory, and the evolution of the luminosity 
function is determined observationally, this connection can 
provide a constraint on the evolution of the LV relation. 

Determining the number density of galaxies as a func- 
tion of their rest-frame luminosity is considerably easier 
than measuring their rotation velocities directly at high 
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redshift. Observed luminosities can be k-corrected to the 
rest visual frame up to z ~ 2 using multi-band photome- 
try extending into the near-IR. It has been very difficult 
to obtain spectroscopic redshifts in the 'desert' regime of 
1.4 ;$ z £j 2, but photometric redshifts offer a means of fill- 
ing in this gap. These advances, combined with theoretical 
progress in predicting the number density and structural 
properties of dark matter halos, suggest the approach for 
constraining the redshift evolution of the LV relation that 
we exploit here. 

The basis of our analysis is the evolution of the velocity 
function of DM halos as determined via a high-resolution 
N-body simulation of the ACDMcosmology (Sigad et al. 
2000, hereafter S00). We extend these results to other 
cosmologies using analytic models. We consider a wide 
range of possibilities for the redshift evolution of the LV 
relation, including some motivated by the results of semi- 
analytic models (cf. Somerville & Primack 1999). We then 
present predictions of the corresponding redshift evolution 
of the comoving number density of bright galaxies, which 
may be compared directly with observations to provide 
quantitative constraints on how the LV relation evolves in 
time. 



2. VELOCITY FUNCTION 

We have calculated the velocity function of DM halos 
using an N-body simulation of the currently popular flat 
ACDM model with ft m = 0.3, fl A = 0.7, Hubble con- 
stant h = 0.7, and ag = 1.0. The simulation (Klypin et 
al. 1999) used a high-resolution adaptive refinement tree 
code (Kravtsov, Klypin & Khokhlov 1997) to follow the 
evolution of structure in a periodic box of comoving side 
60ft. _1 Mpc, with a force resolution of ~ 2h~ 1 kpc in the 
densest regions. 

We identify bound, virialized halos using a spherical 
overdensity-based algorithm and the standard definition 
of virial density from the spherical top-hat collapse model 
(cf. Eke, Cole, & Frenk 1996; Bryan & Norman 1999). 
Details of the halo finder are described in Bullock et 
al. (2000; hereafter BOO). Although the virial velocity 
v v = [67m v /r v ] 1 / 2 (where m v and r v ar the virial mass 
and radius) is often used to approximate galactic circular 
velocities, a more relevant quantity observationally is the 
maximum rotation velocity of the halo v m , which can be 
related to v v via the halo density profile. Over relevant 
radii, halo density profiles are well approximated by the 
NFW functional form with two free parameters (Navarro, 
Frenk & White 1996): p{r) = p^r/r^-^l + r/r s )- 2 , 
where r s is the radius at which the log slope of the pro- 
file is —2. The halo maximum circular velocity occurs at 
2.15r s and is related to v v via v m — \J c//(c)w v , where 
/(c) = 4.62[ln(l + c) - c(l + c)" 1 ]. Here, c is the "con- 
centration" parameter c = r v /r s . The value of c is closely 
related to the halo mass, and for a typical DM halo sur- 
rounding a bright galaxy (m v ~ lO 12 /^ 1 Ai©) at z = 0, 
our simulations yield c ~ 13, implying v m — 1.3u v . The 
concentration at a fixed mass evolves rapidly with redshift, 
c oc (1 + z)~ l (BOO), implying a corresponding evolution 
in v m (v v ). The halo catalogs described in BOO provide v m 
and v v for each halo. A detailed study of the halo velocity 



function derived from these catalogs is presented by S00. 
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Fig. 1. — Evolution of relative comoving number density 
for fixed u m = 200 km s -1 (bold curves) and v v = 200kms _1 
halos in three cosmologies. 

Figure 1 shows the evolution with redshift of the comov- 
ing number density of simulated halos of fixed velocity, ei- 
ther v m = 200 kms -1 or v v — 200 kms -1 . The evolution 
of halos of fixed v m is weaker than for halos of fixed v v due 
to the strong evolution of c. A fit 



dn = V>(^m)d(logv m ) = V>*i>md(logi; m ) 



(1) 



for the whole halo v m distribution is provided by S00, 
who find that the comoving density at v m = 300 km s -1 
is constant as a function of z (±0.1 in the log) at least 
out to z = 5, while (3 w -2.9 - 2.05 + 1.4z 2 ± 0.1, where 
z = log(l + z). 

Figure 1 also shows predictions based on the analytic ap- 
proximation of Sheth & Tormen (1999, hereafter ST) — a 
modification of the Press-Schechter approximation (1974) 
which improves the agreement with simulations. We con- 
verted the ST mass to v v and then to v m using the analytic 
model of BOO for c(m v , z). We see that for ACDM, the an- 
alytic model provides a good approximation to the simu- 
lated evolution both for fixed v m and fixed v v . This allows 
us to extend our predictions using ST to two other cosmo- 
logical models: rCDM, a flat model with Q m = 1, h = 0.5 
and a power spectrum of shape parameter T = 0.21, and 
OCDM, an open geometry with Q ro = 0.3, Ha — 0, and 
h = 0.7. 

In the case of ACDM, we note an increase of ~ 30% in 
the density for constant v m between z — and z ~ 3, and 
a gradual decline back to today's density by z ~ 5. The 
density of halos of a constant v v grows by a factor ~ 2.6 
by z ~ 3, and then declines at higher z. Note that a halo 
with a fixed v v has a decreasing mass and virial radius 
at increasing redshift. Compared to ACDM, the OCDM 
model predicts a slightly larger increase in the density of 
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200 km s" 1 halos at all redshifts. In the case of rCDM, the 
peak density is obtained at z ~ 1, followed by a sharper 
decline, but note that even in this case the comoving den- 
sity at z ~ 2 is comparable to the density at z = 0. 

3. LV RELATIONS AND LUMINOSITY FUNCTIONS 

Armed with this knowledge about the evolution of the 
halo velocity function, we can now predict the evolution 
of the number density of galaxies brighter than a fixed 
absolute luminosity for a given model of LV evolution. 

A possible hitch is that the halo maximum velocity v m 
typically occurs at a radius r m that lies outside the lumi- 
nous parts of the galaxy for which the LV quantities are 
measured. For example, a v m = 200 kms -1 halo at z = 
has r m ~ 45kpc. The velocity due to the halo alone obvi- 
ously declines at smaller radii, but the condensed baryons 
contribute significantly to the total circular velocity there. 
We therefore make the assumption that the rotation curve 
is roughly flat between where it is measured and r m , i.e. 
that the halo v m is the same as the velocity used in the 
observed LV relations. Gonzalez et al. (1999; G99) 5 used 
this assumption at z = to compare the velocity function 
of DM halos with that obtained from observed luminos- 
ity functions and LV relations. The good agreement that 
they obtained is evidence that this is a reasonable approx- 
imation at low redshift. We now assume that this relation 
holds at high redshift as well. 

We parameterize the mean LV relation as: Mb — 
a(z) — b(\og 2v m — 2.5) + 5 log/i, where only the zero-point 
may evolve in time according to the parametric functional 
form 



a(z) =a -2.5 logics), F(z) = (l + z a ) 



0/c 



(2) 



We choose ao = —18.71 and b = 6.76, corresponding to the 
empirical local LV relation used in G99. The parameters 
a and j3 allow a wide range of shapes, where a governs 
the behavior at z < 1, and j3 the asymptotic shape at 
larger redshifts. The indications for only little evolution 
at z £ 1 guide us to consider mostly large positive values 
for a, while the value of (3 is unconstrained, and even its 
sign is unknown. 

In Figure 2 we consider five different representative cases 
for the LV zero-point evolution a{z), and show the pre- 
dicted evolution of the number density of galaxies with 
Mb < —20.2, relative to their number density at z = 0. 
The case (a, 0) — (0, 0) is the case of no evolution in the 
LV relation. The case (2, —1) is close to the evolution pre- 
dicted for the ACDM model if the mass-to-light ratio is 
constant in time. In this model, because halos of fixed 
v m have smaller masses at high redshift, they are fainter 
and therefore the predicted number of bright galaxies de- 
creases. 

The case (1,1) is close to the prediction of a semi- 
analytic model (the "accelerated quiescent" model of 
Somerville, Primack, & Faber 2000) in which the efficiency 
of star formation scales with the inverse dynamical time of 
the halo. This scaling is motivated by one of the empirical 
recipes suggested Kennicutt (1998), based on observations 

5 The published version of this manuscript contains an error in the 



of nearby galaxies. The dynamical time for a halo of a fixed 
velocity is smaller at high redshift (because of the higher 
density), and so galaxies of a fixed velocity are brighter. 
This model predicts an increase in the number density of 
bright galaxies out to z ~ 5. Although this star forma- 
tion recipe is a standard choice in semi-analytic models 
(e.g., Kauffmann, White, & Guiderdoni 1993), note that 
it predicts fairly strong evolution in the LV relation even 
at z <~ 1. The other two cases, (5, 1) and (4, —4), are ar- 
bitrary choices that illustrate a wide range of possibilities 
while keeping the LV relation almost unchanged out to 
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Fig. 2. — Five parameterized models marked by (a, j3) 
for the evolution of the LV zero-point, ao — a(z) (Eq. 1; 
larger values are brighter), along with the implied relative 
evolution of the comoving density of galaxies brighter than 
Mb = —20.2 for the three cosmological models. 



4. DISCUSSION 

We have argued that observed estimates of the number 
density of galaxies above a given luminosity as a function 
of redshift can be used to constrain the redshift evolu- 
tion of the LV relation. Lilly ct al. (1995) found only 
little evolution in the number density of bright galaxies 
in the rest-frame B-band out to z ~ 1, consistent with 
the direct constraints on LV evolution out to this red- 
shift (Vogt et al. 1996, 1997, 2000). Redshift surveys 
of the past decade were unable to reach the higher red- 
shift range 1.0 ^ z ^ 2.0 associated with the spectroscopic 
'desert' or 'high place of sacrifice'. The results presented 
by Dickinson (2000), based on NICMOS (near-IR) imag- 
ing of the Hubble Deep Field (HDF) suggest that the as- 
sembly of present-day luminous galaxies may be taking 
place precisely in this redshift interval. Using photometric 
redshifts, Dickinson finds that there are only one-third as 
many galaxies with rest-frame My < —19 in the redshift 

Equation 12: the quantity in the square root should be inverted. 
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interval 1.4 ^ z 2 as in the interval £ z ;$ 1.4. In our 
fiducial ACDM cosmology, the comoving volumes of the 
HDF for these two intervals are equal, implying an overall 
drop in the comoving number density at z ;> 1.4 of about 
a factor of three. Taken at face value, this would corre- 
spond to dramatic evolution in the LV relation at z J> 1, 
rather similar to the most extreme 'dimming' model shown 
in Figure 2. 

Probably the weakest link in our calculation is our as- 
sumption that the halo maximum velocity is the same as 
the measured internal velocity for a galaxy. The actual 
measured velocity will depend on the baryon fraction in 
the galaxy, the spatial distribution of the baryons, and 
the degree to which angular momentum is preserved in 
the collapse. However, note that because we expressed 
the space density in Figure 2 relative to that at z = 0, 
our plotted results will not change by much as long as the 
ratio of observed velocity to the halo maximum velocity, 
C = u bs/"m, is of order unity and roughly constant with 
redshift. We have calculated the value of C(z) using a 
simple analytic model of disk formation, which assumes 
that a constant fraction of the halo mass condenses into 
a disk, and angular momentum is conserved (Mo, Mao,& 
White 1998). We find that C{z) ~ 1.3 for the redshift 
range < z < 2, and falls to C{z) ~ 1.0 by z ~ 5. This 
suggests that our predictions should be reasonably robust 
at least out to z ~ 2. 

How significant is the effect seen in the HDF? One pos- 
sibility is that galaxies are preferentially missed at z ;> 1.4 
because of the (l+z) 4 cosmological surface-brightness dim- 
ming. Dickinson (2000) claims that this is unlikely to be 
a large effect; when bright HDF galaxies from z ;$ 1 are 
artificially rcdshifted to z ^ 2 with no intrinsic luminosity 
evolution, most of them are still easily detectable. An- 
other possibility is that the photometric redshifts are in- 
accurate; Dickinson points out that there are no spectro- 
scopic redshifts in the range 1.4 z ,$ 2 where the deficit 
is seen. However, this would require quite a dramatic fail- 



ure where all of the galaxies (with photometric redshifts) 
assigned to 0.5 z p hot ~ 1-4 would have to be reassigned 
to the higher redshift range. Most importantly, recall that 
the HDF covers a relatively small volume, and therefore 
random fluctuations along the line of sight are likely to be 
large. The observed deficit of galaxies might simply reflect 
large-scale "spikes" in the redshift distribution. Although 
there is some evidence for a similar deficit of bright galax- 
ies at z £ 1 in other deep optical/IR surveys (Fontana et 
al. 1999), secure confirmation of these results will require 
deep wide-field surveys with multi-band optical-IR pho- 
tometry. Several such surveys are now in progress (e.g. the 
Cambridge-Carnegie Las Campanas IR Survey). More- 
over, the DEEP survey (Davis & Faber 1998) will obtain 
measured line-widths for a large number of galaxies at 
z J> 0.7, which will provide direct constraints on the LV 
relation at high redshift. 

Finally, our results may be relevant to a recent proposal 
by Newman & Davis (2000) to use velocity-selected galax- 
ies as a tracer population in a classical alN/dz measure of 
cosmological parameters. This analysis would be subject 
to similar uncertainties in relating the halo velocity to ob- 
servable galaxy velocity, e.g., the dependence on whether 
the halos are selected by v m or v v . Moreover, because the 
proposed spectroscopic surveys are magnitude limited, the 
indicated evolution of the LV relation could introduce se- 
vere selection biases into the analysis. 
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